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Nonequilibrium. 



Systems in stationary nonequilibrium are mechanical systems subject to nonconservative ex- 
ternal forces and to thermostat forces which forbid indefinite increase of the energy and allow 
reaching statistically stationary states. A system £ is described by the positions and velocities of 
its n particles X, X, with the particles positions confined to a finite volume container Co- 

If X = (xi, . . ., x n ) are the particles positions in a Cartesian inertial system of coordinates, 
the equations of motion are determined by their masses m; > 0, i = l,...,n, by the potential 
energy of interaction V(xi, . . . ,x n ) = V(X.), by the external nonconservative forces Fi(X,^) and 
by the thermostat forces — i?j as 

m i Z i = -d Xi V(X) + F i {X;&)-'& i , i = l,...,n (1.1) 

where # = {fx, ■ ■ ■ , f q ) are strength parameters on which the external forces depend. All forces and 
potentials will be supposed smooth, i.e. analytic, in their variables aside from possible impulsive 
elastic forces describing shocks, and with the property: F(X; 0) = 0. The impulsive forces are 
allowed here to model possible shocks with the walls of the container Co or between hard core 
particles. 

A kind of thermostats are reservoirs which may consist of one or more infinite systems which 
are asymptotically in thermal equilibrium and are separated by boundary surfaces from each other 
aas well as from the system: with the latter they interact through short range conservative forces, 
see Fig.l. 

The reservoirs occupy infinite regions of the space out- 




1 3 

side Co, e.g. sectors C a C R , a = 1, 2 . . ., in space and 
their particles are in a configuration which is typical of 
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an equilibrium state at temperature T a . This means that 
X Fig 1 ^ e empirical probability of configurations in each C a is 

Gibbsian with some temperature T a . In other words the 
frequency with which a configuration (Y, Y + r) occurs in 
a region A + r C C a and a configuration (W, W + r) occurs outside A + r (with Y C A, W n A 
= 0) averaged over the translations A + r of A by r (with the restriction that A + r C C a ) is 

e -/3a(d^|Y| 2 +V Q (Y|W)) 

average (/ A+r [(Y, Y + r); W, W + r]) = — (1.2) 

r _l_AcC a normalization 
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here m a is the mass of the particles in the a-th reservoir and V a (Y|W) is the energy of the short 
range potential between pairs of particles in Y C C a or with one point in Y and one in W. 
Since the configurations in the system and in the thermostats are not random the (1.2) should be 
considered as an "empirical" probability in the sense that it is the frequency density of the events 
{(Y, Y + r); W + r}: in other words the configurations u a in the reservoirs should be "typical" in 
the sense of probability theory of distributions which are asymptotically Gibbsian. 

The property of being "thermostats" means that (1.2) remains true for all times, if initially 
satisfied. 

Mathematically there is a problem at this point: the latter property is either true or false, 
but a proof of its validity seems out of reach of the present techniques except in very simple 
cases. Therefore here we follow an intuitive approach and assume that such thermostats exist and, 
actually, that any configuration which is typical of a stationary state of an infinite size system of 
interacting particles in the C 's, with physically reasonable microscopic interactions, satisfies the 
property (1.2). 

The above thermostats are examples of "deterministic thermostats" because, together with 
the system they form a deterministic dynamical system. They are called "Hamiltonian thermostats" 
and are often considered as the most appropriate models of "physical thermostats" . 

A closely related thermostat model is obtained by assuming that the particles outside the 
system are not in a given configuration but they have a probability distribution whose conditional 
distributions satisfy (1.2) initially. Also in this case it is necessary to assume that (1.2) remains 
true for all times, if initially satisfied. Such thermostats are examples of "stochastic thermostats" 
because their action on the system depends on random variables u a which are the initial configu- 
rations of the particles belonging to the thermostats. 

Other kinds of stochastic thermostats are collision rules with the container boundary dCo 
of S: every time a particle collides with dC® it is reflected with a momentum p in d 3 p that 
has a probability distribution proportional to e~@ a ^ p d 3 p where a , a = 1, 2, . . . depend on which 
boundary portion (labeled by a = 1, 2 . . . and, if ks is Boltzmann's constant, at "temperature" T a = 
(fcs/Sa) -1 ) the collision takes place. Which p is actually chosen after each collision is determined 
by a random variable u> = (u>i, u> 2 , ...). 

It is also possible, and convenient, to consider deterministic thermostats which are finite. In 
the latter case # is a force only depending upon the configuration of the n particles in their finite 
container C . The distinction between stochastic and deterministic thermostats ultimately rests on 
what we call "system" . If reservoirs or the randomness generators are included in the system then 
the system becomes deterministic (possibly infinite); and finite deterministic thermostats can also 
regarded as simplified models for infinite reservoirs, see Sect. 4. 

Examples of finite deterministic reservoirs are forces obtained by imposing a nonholonomic 
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constraint via some ad hoc principle like the Gauss' principle. For instance if a system of particles 

def 

driven by a force Gi = —d Xi V(X.) + Fj (X) is enclosed in a box Co and # is a thermostat enforcing 
an anholonomic constraint Y>(X, X) = via Gauss' principle then 



0i(x,x)= — 3 \ ' m 3 — 3 — 1 a x .v(x,x) (i.3) 

Gauss' principle is remarkable as it says that the force which needs to be added to the other forces 
Gi acting on the system minimizes J2i ( Gi ~™ i3 -^ , given X, X, among all accelerations which 
are compatible with the constraint ip. 

For simplicity stochastic or infinite thermostats will not be considered here. It should be kept 
in mind that the only known examples of mathematically treatable thermostats modeled by infinite 
reservoirs are cases in which the thermostats particles are either noninteracting particles or linear 
(i.e. noninteracting) oscillators. 

In general in order that a force i9 can be considered a deterministic "thermostat force" a further 
property is necessary: namely that the system evolves according to (1.1) towards a stationary 
state. This means that for all initial particles configurations (X, X), except possibly for a set of 
zero phase space volume, any smooth function /(X, X) evolves in time so that, if S t (X, X) denotes 
the configuration into which the initial data evolve in time t according to (1.1), then the limit 

T 

exists and is independent of (X, X). The probability distribution /x is then called the SRB distri- 
bution for the system. The maps S t will have the group property S t ■ S t > = S t+t > and the SRB 
distribution fj, will be invariant under time evolution. 

It is important to stress that the requirement that the exceptional configurations form just a 
set of zero phase volume (rather than a set of zero probability with respect to another distribution, 
singular with respect to the phase volume) is a strong assumption and it should be considered an 
axiom of the theory: it corresponds to the assumption that the initial configuration is prepared 
as a typical configuration of an equilibrium state, which by the classical equidistribution axiom of 
equilibrium statistical mechanics is a typical configuration with respect to the phase volume . 

For this reason the SRB distribution is said to describe a stationary state of the system. The 
SRB distribution depends on the parameters on which the forces acting on the system depend, 
e.g. | Co | (volume), ^ (strength of the forcings), (temperatures) etc. The collection of SRB 

distributions obtained by letting the parameters vary defines a nonequilibrium ensemble. 

In the stochastic case the distribution \i is required to be invariant in the sense that it can 
be regarded as a marginal distribution of an invariant distribution for the larger (deterministic) 
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system formed by the thermostats and the system itself. 
References: [EM90], [Ru97a], [EPR99]. 

2. Nonequilibrium thermodynamics 

The key problem of nonequilibrium statistical mechanics is to derive a macroscopic "nonequi- 
librium thermodynamics" in a way similar to the derivation of equilibrium thermodynamics from 
equilibrium statistical mechanics. 

The first difficulty is that nonequilibrium thermodynamics is not well understood. For instance 
there is no (agreed upon) definition of entropy, while it should be kept in mind that the effort to 
find the microscopic interpretation of equilibrium entropy, as defined by Clausius, was a driving 
factor in the foundations of equilibrium statistical mechanics. 

The importance of entropy in classical equilibrium thermodynamics rests on the implication 
of universal, parameter free, relations which follow from its existence (e.g. dy^ = du^ if U is the 
internal energy, T the absolute temperature and p the pressure of a simple homogeneous material). 

Are there universal relations among averages of observables with respect to SRB distributions? 

The question has to be posed for systems "really" out of equilibrium, i.e. for ^ 7^ (see 
(1.1)): in fact there is a well developed theory of the derivatives with respect to ^ of averages of 
observables evaluated at ^ = 0. The latter theory is often called, and here we shall do so as well, 
classical nonequilibrium thermodynamics or near equilibrium thermodynamics and it has been quite 
successfully developed on the basis of the notions of equilibrium thermodynamics, paying particular 
attention on the macroscopic evolution of systems described by macroscopic continuum equations 
of motion. 

"Stationary nonequilibrium statistical mechanics" will indicate a theory of the relations be- 
tween averages of observables with respect to SRB distributions. Systems so large that their 
volume elements can be regarded as being in locally stationary nonequilibrium states could also be 
considered. This would extend the familiar "local equilibrium states" of classical nonequilibrium 
thermodynamics: however they are not considered here. This means that we shall not attempt at 
finding the macroscopic equations regulating the time evolution of continua locally in nonequilib- 
rium stationary states but we shall only try to determine the properties of their "volume elements" 
assuming that the time scale for the evolution of large assemblies of volume elements is slow com- 
pared to the time scales necessary to reach local stationarity. 

References: [DGM84], [Le93], [Ru97a], [Ru99a], [Ga98], [GL03], [Ga04]. 

3. Chaotic hypothesis 

In equilibrium statistical mechanics the ergodic hypothesis plays an important conceptual role 
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as it implies that the motions of ergodic systems have a SRB statistics and that the latter coincides 
with the Liouville distribution on the energy surface. 

An analogous role has been proposed for the chaotic hypothesis: which states that the motion 
of a chaotic system, developing on its attracting set, can be regarded as an Anosov system. This 
means that the attracting sets of chaotic systems, physically defined as systems with at least 
one positive Lyapunov exponent, can be regarded as smooth surfaces on which motion is highly 
unstable: 

(i) around every point a curvilinear coordinate system can be established which has three planes, 
varying continuously with x, which are covariant (i.e. are coordinate planes at a point x which are 
mapped, by the evolution St, into the corresponding coordinate planes around Stx) and 

(ii) the planes are of three types, stable, unstable and marginal, with respective positive dimensions 
d s ,d u and 1: lengths on the stable surface and on the unstable surface of any point contract at 
exponential rate as time proceeds towards the future or towards the past. The length along the 
marginal direction neither contracts nor expands (i.e. it varies around the initial value staying 
bounded away from and oo): its tangent vector is parallel to the flow. In cases in which time 
evolution is discrete, and determined by a map S, the marginal direction is missing. 

(iii) the contraction over a time t, positive for lines on the stable plane and negative for those on 
the unstable plane, is exponential, i.e. lengths are contracted by a factor uniformly bounded by 
Ce- K l*l with C,k > 0. 

(iv) there is a dense trajectory. 

It has to be stressed that the chaotic hypothesis concerns physical systems: mathematically 
it is very easy to find dynamical systems for which it does not hold. As it is easy (actually even 
easier) to find systems in which the ergodic hypothesis does not hold (e.g. harmonic lattices or 
black body radiation). However, if suitably interpreted, the ergodic hypothesis leads even for these 
systems to physically correct results (the specific heats at high temperature, the Raileigh-Jeans 
distribution at low frequencies). Moreover the failures of the ergodic hypothesis in physically 
important systems have led to new scientific paradigms (like quantum mechanics from the specific 
heats at low temperature and Planck's law). 

Since physical systems are almost always not Anosov systems it is very likely that probing 
motions in extreme regimes will make visible the features that distinguish Anosov systems from 
non Anosov systems: much as it happens with the ergodic hypothesis. 

The interest of the hypothesis is to provide a framework in which properties like the existence 
of an SRB distribution is a priori guaranteed: the role of Anosov systems in chaotic dynamics is 
similar to the role of harmonic oscillators in the theory of regular motions. They are the paradigm 
of chaotic systems as the harmonic oscillators are the paradigm of order. Of course the hypothesis 



2/ febbraio/2008; 17:31 



5 



is only a beginning and one has to learn how to extract information from it, as it was the case 
with the use of the Liouville distribution once the ergodic hypothesis guaranteed that it was the 
appropriate distribution for the study of the statistics of motions in equilibrium situations. 

References: [Ru76], [GC95], [Ru97a], [Ga98], [GBG04]. 

4. Heat, temperature and entropy production 

The amount of heat Q that a system produces while in a stationary state is naturally identified 
with the work that the thermostat forces # perform per unit time 

i 

A system may be in contact with several reservoirs: in models this will be reflected by a decompo- 
sition 

# = ^tf( a )(X,X) (4.2) 

a=l 

where is the force due to the a-th thermostat and depends on the coordinates of the particles 
which are in a region A a C Co of a decomposition \J™ =1 A a = Co of the container Co occupied by 
the system (A a n A a / = if a ^ a'). 

From several studies based on simulations of thermostated systems of particles arose the pro- 
posal to consider the average of the phase space contraction cr( a )(X, X) due to the a-th thermostat 



a(a) ^ x) d g ^ 0( O) ( ^ x) (43) 



3 

and to identify it with the rate of entropy creation in the a-th thermostat. 

Another key notion in thermodynamics is the temperature of a reservoir; in the infinite de- 
terministic thermostats case, of Sect.l, it is defined as (ksPa) -1 but in the finite determinstic 
thermostats considered here it needs to be defined. If there are m reservoirs with which the system 
is in contact one sets 

a (a) dg ^ (a) ^ x ^ ^ j a(a) ^ x) ^ (dX dX) 

q.^y:*p-** (4 ' 4) 

i 

where \i is the SRB distribution describing the stationary state. It is natural to define, if ks is 
Boltzmann's constant, the absolute temperature of the a-th thermostat to be 
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T a = ™ (4.5) 

Although it is known that ^ a > it is not clear that T a > 0: this happens in a rather 
general class of models and it would be desirable, for the interpretation that is proposed here, that 
it could be considered a property to be added to the requirements that the forces i? a be thermostats 
models. 

An important class of thermostats for which the property T a > holds can be described as 
follows. Imagine N particles in a container Co interacting via a potential Vo = £i<j ¥>(<li — + 
^"'(Qj) (where V models external conservative forces like obstacles, walls, gravity, . . .) and, 
furthermore, interacting with M other systems S a , of N a particles of mass m a , in containers C a 
contiguous to Co- The latter will model M parts of the system in contact with thermostats at 
temperatures T a , a = 1, . . . , M. 

The coordinates of the particles in the a-th system S a will be denoted x", j = 1, . . . , N a , and 
they will interact with each other via a potential V a = Ylfj <^a(x" — x"). Furthermore there will 
be an interaction between the particles of each thermostat and those of the system via potentials 
W a = Eti £f=i - x|), a = 1, . . . , M. 

The potentials will be assumed to be either hard core or non singular potentials and the 
external V is supposed to be at least such that it forbids existence of obvious constants of motion. 

The temperature of each S a will be defined by the total kinetic energy of its particles, i.e. by 
K a = Ylf=i \ m a{^- a j) 2 d — \N a k,BT a : the particles of the a-th thermostat will be kept at constant 
temperature by further forces t?|. The latter are defined by imposing via a Gaussian constraint 
that K a is a constant of motion (see (1.3) with i/j = K a ) . This means that the equations of motion 
are 



™* = (MQ) + E W a(Q^ a )) 

a=l ( 4 - 6 ) 

m a ^ = -d x a (F a (x a ) + W (Q, x a )) - ^ 

and an application of Gauss' principle yields i9| = 3 ^ Q ^^t "= where L a is the work per 

unit time done by the particles in Co on the particles of S a and V a is their potential energy. 

In this case the partial divergence a a = 3N a a a = k^r~ wm ma ke (4.5) identically 

satisfied with T a > because L a can be naturally interpreted as heat Q a ceded, per unit time, by 
the particles in Co to the subsystem S a (hence to the a-th thermostat because the temperature of 
S a is constant), while the derivative of V a will not contribute to the value o/cr" . The phase space 
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contraction rate is, neglecting the total derivative terms, 

N a ^ 

a true (X,X) = ^-^. (4.7) 

a=l 

where the subscript "true" is to remind that an additive total derivative term distinguishes it from 
the complete phase space contraction. 

Remarks: (1) The above formula provides the motivation of the name "entropy creation rate" 
attributed to the phase space contraction a. Note that in this way the definition of entropy creation 
is "reduced" to the equilibrium notion because what is being defined is the entropy increase of the 
thermostats which have to be considered in equilibrium. No attempt is made here to define the 
entropy of the stationary state. Nor any attempt is made to define the notion of temperature 
of the nonequilibrium system in Cq (the T a are temperatures of the S a , not of the particles in 
Co). This is an important point as it leaves open the possibility of envisaging the notion of "local 
equilibrium" which becomes necessary in the approximation (not considered here) in which the 
system is regarded continuum. 

(2) In the above model another viewpoint is possible: i.e. to consider the system to consist of only 
the N particles in Co and the M systems S a to be thermostats. From this point of view the above 
can be considered a model of a system subject to thermostats. The Gibbs distribution characterizing 
the infinite thermostats of Sect.l becomes in this case the constraint that the kinetic energies K a 
are constants, enforced by the Gaussian forces. This shows that the phase space contraction can be 
an appropriate definition of entropy creation rate only if the system is subject to finite deterministic 
thermostats. In the new viewpoint the appropriate definition should be simply the r.h.s. of (4.7), 
i.e. the work per unit time done by the forces of the system on the thermostats divided by the 
temperature of the thermostats. This suggests a more general definition of entropy creation rate, 
applying also to thermostats that are often considered "more physical" and that needs to be further 
investigated. 

References: [EM90], [GC95], [Ru96], [Ru97b], [Ga04]. 

5. Thermodynamic fluxes and forces 

Nonequilibrium stationary states depend upon external parameters ipj like the temperatures 
T a of the thermostats or the size of the force parameters ^ = ((f\, . . . ,<p q ), see (1.1). Nonequi- 
librium thermodynamics is well developed at "low forcing": strictly speaking this means that it 
is widely believed that we understand properties of the derivatives of the averages of observables 
with respect to the external parameters if evaluated at ipj = 0. Important notions are the notions 
of thermodynamic fluxes Ji and of thermodynamic forces ipi; hence it seems important to extend 
such notions to nonequilibrium systems (i.e. & ^ 0). 
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A possible extension could be to define the thermodynamic flux Jj associated with a force ipi 
as Jj = (d Vi a) SRB where cr(X, X; is the volume contraction per unit time. This definition seems 
appropriate in several concrete cases that have been studied and it is appealing for its generality. 

An interesting example is provided by the model of thermostated system in (4.6): if the 
container of the system is a box with periodic boundary conditions one can imagine to add an 
extra constant force E acting on the particles in the container. Imagining the particles to be 
charged by a charge e and regarding such force as an electric field the first equation in (4.6) is 
modified by the addition of a term eE. 

The constraints on the thermostats temperatures imply that a depends also on E: in fact, if 
J = eJ2j is the electric current, energy balance implies U to t = E • J — ^ a {L a — V a ) if U to t is 
the sum of all kinetic and potential energies. Then the phase space contraction ^ a La ~ Va can be 
written, to first order in the temperature variations 5T a with respect to a common value T a = T, 
as - E a L ^ L§ ^ L + E,J r C/tot hence a true , see (4.7), is 

E J Q a 5T a 

a 

The definition and extension of the conjugacy between thermodynamic forces and fluxes is 
compatible with the key results of classical nonequilibrium thermodynamics, at least as far as 
Onsager reciprocity and Green-Kubo's formulae are concerned. It can be checked that if the 
equilibrium system is reversible, i.e. if there is an isometry / on phase space which anticommutes 
with the evolution (IS t = S- t I in the case of continuous time dynamics t — > St or IS = S _1 I in 
the case of discrete time dynamics S) then, shortening (X, X) into x, 



def 

Lij = d&.Jj\&=o = d^ i (d^ j a(x;^)} SRB \&=o = d*jJi |#=o = 

(5-2) 

dt 

#=o 



1 f°° 

L H=o (d<s>.a{S t x-&)d< s>i (j{x-&)) 

* J -oo 



SRB 



The a(x; <P) plays the role of "Lagrangian" generating the duality between forces and fluxes. The 
extension of the duality just considered might be of interest in situations in which # / 0. 

References: [DGM84] , [Ga96] , [GR97] . 

6. Fluctuations 

As in equilibrium, large statistical fluctuations of observables are of great interest and al- 
ready there is, at the moment, a rather large set of experiments dedicated to the analysis of large 
fluctuations in stationary states out of equilibrium. 

If one defines the dimensionless phase space contraction 
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„(*) = ! f^Mrft (6.1) 
T Jo °+ 

(see also (4.7)) then there exists p* > 1 such that the probability P T of the event p 6 [a, 6] with 
[a, 6] C (— p*,p*) has the form 

P T (j)£ [a, 6]) = const e ™ax p6f0)i)] C( P )+0(i) ( 6 _ 2 ) 

with C(p) analytic in (— p*,p*). The function £(p) can be conveniently normalized to have value 
at p = 1 (i.e. at the average value of p). 

Then, m Anosov systems which are reversible and dissipative (see Sect. 5) a general symme- 
try property, called the fluctuation theorem and reflecting the reversibility symmetry, yields the 
parameterless relation 



C(-P) = C(P)-P*+ (6-3) 

This relation is interesting because it has no free parameters, in other words it is universal for 
reversible dissipative Anosov systems. In connection with the duality fluxes-forces in Sec. 5, it can 
be checked to reduce to the Green-Kubo formula and to Onsager reciprocity, see (5.2), in the case 
in which the evolution depends on several fields ^ and ^ — > (of course the relation becomes 
trivial as # — > because a+ — > and to obtain the result one has first to divide both sides by 
suitable powers of the fields 

A more informal (but imprecise) way of writing (6. 2), (6. 3) is 

Pt ^ - e Tpa + +0 ^\ for all p G (- P *,p*) (6.4) 



PA-P) 

where P T (p) is the probability density of p. An interesting consequence of (6.4) is (c~ T pa+ ) SRB = 1 
in the sense that ^ log (e~ Tpcr+ } SRB T ^oo> 0. 

Occasionally systems with singularities have to be considered: in such cases the relation (6.3) 
may change in the sense that the function ((p) may be not analytic: in such cases one expects that 
the relation holds in the largest analyticity interval symmetric around the origin. In various cases 
considered in the literature such interval appears to contain the interval (—1,1) and sometimes this 
can be proved rigorously. For instance in simple, although admittedly special, examples of systems 
close to equilibrium. 

It is important to remark that the above fluctuation relation is the first representative of 
remarkable consequences of the reversibility and chaotic hypotheses. For instance given Fi, . . . ,F n 
arbitrary observables which are (say) odd under time reversal I (i.e. F(Ix) = —F(x)) and given 
n functions t G [— ^, |] — > <fj(t), j = 1, . . . ,n one can ask which is the probability that Fj(S t x) 
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"closely follows" the pattern tpj (t) and at the same time \ Jq dO has value p. Then calling 

P T {F\ ~ tpi, . . . ,F n ~ (p n ,p) the probability of this event, which we write in the imprecise form 

def 

corresponding to (6.4) for simplicity, and defining Iifj(t) = — ifj(-t) it is 

P T (iWi,---,^~^,p) _„ T<r+ „ pe{ _ p * }pn (6 .5) 



P T (Fi ~ Iip 1: . . . , F n ~ Itf n , -p) 
which is remarkable because it is parameterless and at the same time surprisingly independent of 
the choice of the observables Fj. The relation (6.5) has far reaching consequences: for instance 
if n = 1 and F\ = <9<j,.ct(x; #) the relation (6.5) has been used to derive the mentioned Onsager 
reciprocity and Green-Kubo's formulae at # = 0. 

Eq. (6.5) can be read as follows: the probability that the observables Fj follow given evolution 
patterns (fj conditioned to entropy creation rate pcr+ is the same that they follow the time reversed 
patterns if conditioned to entropy creation rate — pa+ . In other words to change the sign of time 
it is just sufficient to reverse the sign of entropy creation rate, no "extra effort" is needed. 

References: [Si72] , [Si94] , [ECM93] , [GC85] , [Ga96] , [GR97] , [Ga99] , [GBG04] , [BGGZ05] . 

7. Fractal attractors. Pairing. Time reversal. 

Attracting sets (i.e. sets which are the closure of attractors) are fractal in most dissipative 
systems. However the chaotic hypothesis assumes that fractality can be neglected. Aside from the 
very interesting cases in which systems are close to equilibrium, in which the closure of an attractor 
is the whole phase space (under the chaotic hypothesis, i.e. if the system is Anosov) there are, 
however, serious problems in preserving validity of the fluctuation theorem. 

The reason is very simple: if the attractor closure is smaller than phase space then it is to be 
expected that time reversal will change the attractor into a repeller disjoint from it. Thus even if 
the chaotic hypothesis is assumed, so that the attracting set A can be considered a smooth surface, 
the motion on the attractor will not be time reversal symmetric (as its time reversal image will 
develop on the repeller): one can say that an attracting set with dimension lower than that of 
phase space in a time reversible system corresponds to a spontaneous breakdown of time reversal 
symmetry. 

It has been noted however that there are classes of systems, forming a large set in the space of 
evolutions depending on a parameter $, in which geometric reasons imply that if beyond a critical 
value $ c the attracting set becomes smaller than phase space then a map Ip is generated mapping 
the attractor A into the repeller 1Z, and viceversa, such that Ip is the identity on A U 1Z and Ip 
commutes with the evolution: therefore the composition / ■ Ip is a time reversal symmetry {i.e. it 
anticommutes with evolution) for the motions on the attracting set A (as well as on the repeller 
K). 
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In other words the time reversal symmetry in such systems "cannot be broken" : if spontaneous 
breakdown occurs (i.e. A is not mapped into itself under time reversal I) a new symmetry Ip is 
spawned and / ■ Ip is a new time reversal symmetry (an analogy with the spontaneous violation of 
time reversal in quantum theory, where time reversal T is violated but TCP is still a symmetry: 
so T plays the role of I and CP that of Ip). 

Thus a fluctuation relation will hold for the phase space contraction of the motions taking 
place on the attracting set for the class of systems with the geometric property mentioned above 
(technically the latter is called axiom C property). This is interesting but it still is quite far from 
being checkable even in numerical experiments. 

There are nevertheless systems in which also a pairing property holds: this means that, 
considering the case of discrete time maps S, the Jacobian matrix d x S(x) has 2N eigenvalues that 
can be labeled, in decreasing order, Xn(x), . . . , Xi N (x), ...,A\(x) with the remarkable property 
that i(AAr_j(x) + Xj(x)) d = a(x) is j-independent. In such systems a relation can be established 
between phase space contractions in the full phase space and on the surface of the attracting set: the 
fluctuation theorem for the motion on the attracting set can therefore be related to the properties of 
the fluctuations of the total phase space contraction measured on the attracting set (which includes 
the contraction transversal to the attracting set) and if 2M is the attracting set dimension and 2N 
is the total dimension of phase space it is, in the analyticity interval (— p*,p*) of the function ((p), 

C(-p) = C(p)-p^+ (7-1) 

which is an interesting relation. It is however very difficult to test in mechanical systems because 
in such systems it seems very difficult to make the field so high to see an attracting set thinner 
than the whole phase space and still observe large fluctuations. 

References: [DM96], [Ga99]. 

8. Nonequilibrium ensembles and their equivalence 

Given a chaotic system the collection of the SRB distributions associated with the various con- 
trol parameters (volume, density, external forces, ...) forms an "ensemble" describing the possible 
stationary states of the system and their statistical properties. 

As in equilibrium one can imagine that the system can be described equivalently in several 
ways at least when the system is large ("in the thermodynamic" or "macroscopic limit"). In 
nonequilibrium equivalence can be quite different and more structured than in equilibrium because 
one can imagine to change not only the control parameters but also the thermostatting mechanism. 

It is intuitive that a system may behave in the same way under the influence of different 
thermostats: the important phenomenon being the extraction of heat and not the way in which it is 
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extracted from the system. Therefore one should ask when two systems are "physically equivalent" , 
i.e. when the SRB distributions associated with them give the same statistical properties for the 
same observables, at least for the very few observables which are macroscopically relevant. The 
latter may be a few more than the usual ones in equilibrium (temperature, pressure, density, ...) 
and include currents, conducibilities, viscosities, ... but they will always be very few compared to 
the (infinite) number of functions on phase space. 

As an example consider a system of N interacting particles (say hard spheres) of mass m 
moving in a periodic box Co of side L containing a regular array of spherical scatterers (a basic 
model for electrons in a crystal) which reflect particles, elastically and are arranged so that no 
straight line exixts in Co which avoids the obstacles (to eliminate obvious constants of motion) . An 
external field Eu acts also along the u-direction: hence the equations of motion are 



mxi = fi + Eu - (8.1) 

where are the interparticle and scatterers-particles forces and 1?^ are the thermostating forces. 
The following thermostats models have been considered 

(1) I &i = fki {viscosity thermostat) 

(2) immediately after elastic collision with an obstacle the velocity is rescaled to a prefixed value 
yjZkBTm- 1 for some T (Drude's thermostat) 

(3) ■di = r- ^ . 2 ' (Gauss' thermostat) 

The first two are not reversible. At least not manifestly such, because the natural time 
reversal, i.e. change of velocity sign, is not a symmetry (there might be however more hidden, 
hitherto unknown, symmetries which anticommute with time evolution). The third is reversible 
and time reversal is just the change of the velocity sign. The third thermostat model generates a 
time evolution in which the total kinetic energy K is constant. 

Let fi'„, h't, [J>'k be the SRB distributions for the system in a container Co with volume | C 1 = 
L 3 and density p = fixed. Imagine to tune the values of the control parameters v, T, K in such a 
way that (kinetic energy)^ = 8, with the same £ for p = p' u , p'^, p'x and consider a local observable 
F(X, X) > depending only on the coordinates of the particles located in a region A C Co- Then 
a reasonable conjecture is that 

(F) , IF) , 

lim j-p^ = lim -—^ = 1 (8.2) 

N=p K >Ht Mt= P V Vt 

if the limits are taken at fixed F (hence at fixed A while L oo). The conjecture is an 

open problem: it illustrates, however, the a kind of questions arising in nonequilibrium statistical 
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mechanics. 



References: [ES93] , [Ga99] , [Ru99b] . 
9. Outlook 

The subject is (clearly) at a very early stage of development. 

(1) The theory can be extended to stochastic thermostats quite satisfactorily, at least as far 
as the fluctuation theorem is concerned. 

(2) Remarkable works have appeared on the theory of systems which are purely Hamiltonian 
and (therefore) with thermostats that are infinite: unfortunately the infinite thermostats can be 
treated, so far, only if the systems are "free" at infinity (either free gases or harmonic lattices). 

(3) The notion of entropy turns out to be extremely difficult to extend to stationary states 
and there are even doubts that it could be actually extended. Conceptually this is certainly a major 
open problem. 

(4) The statistical properties of stationary states out of equilibrium are still quite mysterious 
and surprising: recently appeared exactly solvable models as well as attempts at unveiling the deep 
reasons for their solubility and at deriving from them general guiding principles. 

(5) Numerical simulations have given a strong impulse to the subject, in fact one can even say 
that they created it: introducing the notion of thermostat and providing the first reliable results 
on the properties of systems out of equilibrium. Simulations continue to be an essential part of the 
effort of research on the field. 

(6) Approach to stationarity leads to many important questions: is there a Lyapunov function 
measuring the distance between an evolving state and the stationary state towards which it evolves? 
In other words can one define an analogous of Boltzmann's i^-function? About this question there 
have been proposals and the answer seems affirmative, but it does not seem that it is possible to find 
a universal, system independent, such function (search for it is related to the problem of defining 
an entropy function for stationary states: its existence is at least controversial, see Sec. 2, 3). 

(7) How irreversible is a given irreversible process in which the initial state fj,o is a stationary 
state and at time t = and the external parameters start changing into functions &(t) of t 
and tend to a limit as t — > oo? In this case the stationary distribution fi starts changing and 
becomes a function fi t of t which is not stationary but approaches another stationary distribution 
/ioo as t — > oo. The process is, in general, irreversible and the question is how to measure its degree 
of irreversibility. A natural quantity I associated with the evolution from an initial stationary state 
to a final stationary state through a change in the control parameters can be defined as follows. 
Consider the distribution jj, t into which /x evolves in time t, and consider also the SRB distribution 
/•*#(*) corresponding to the control parameters "frozen" at the value at time t, i.e. &(t). Let the 
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phase space contraction, when the forces are "frozen" at the value &(t), be a t (x) = a(x;&(t)). In 
general fj, t / H&{t) ■ Then 

T({*(t)}, ^0,^00) = / {Ht{<Tt) ~ H*{t){<?t)) 2 dt (9.1) 
J 

can be called the degree of irreversibility of the process: it has the property that in the limit of 
infinitely slow evolution of &(t), e.g. if &{t) = ^0 + (1 — e _7K ')/\ (a quasi static evolution on time 
scale 7 _1 k _1 from #0 to = ^0 + "^)> the irreversibility degree I 1 Q > if (as in the case 
of Anosov evolutions, hence under the chaotic hypothesis) the approach to a stationary state is 
exponentially fast at fixed external forces ^. 

The entire subject is dominated by the initial insights of Onsager on classical nonequilibrium 
thermodynamics: which concern the properties of the infinitesimal deviations from equilibrium 
(i.e. averages of observables differentiated with respect to the control parameters ^ and evaluated 
at ^ = 0). The present efforts are devoted to studying properties at ^ 7^ 0. In this direction 
the classical theory provides certainly firm constraints (like Onsager reciprocity or Green-Kubo 
relations or fluctuation dissipation theorem) but at a technical level it gives little help to enter the 
terra incognita of nonequilibrium thermodynamics of stationary states. 

References: [Ku98],[LS99],[Ma99]; [EPR99],[BLR00],[EY05]; [DLS01];[BDGJL01]; 
[EM90] , [ECM33] ; [GL03] , [Ga04] , [BGGZ05] . 

Bibliography 

[BDGJL01] L. Bertini, A. De Sole, D. Gabrielli, G. Jona, C. Landim: Fluctuations in Station- 
ary Nonequilibrium States of Irreversible Processes, Physical Review Letters, ??, , 2001.87 (2001) 
040601 

[BLR00] F. Bonetto, J.L. Lebowitz, L. Rey-Bellet: Fourier's law: a challenge to theorists, in 
"Mathematical Physics 2000", p. 128-150, editor R. Streater, Imperial College Press, 2000. 

[BGGZ05] F. Bonetto, G. Gallavotti, A. Giuliani, F. Zamponi: Chaotic Hypothesis, Fluctuation 
Theorem and Singularities, mp_arc 05-257, |cond-mat/0507672| 

[DM96] C. Dettman, G.P. Morriss: Proof of conjugate pairing for an isokinetic thermostat, Physical 
Review 53 E, 5545-5549, 1996. 

[DGM84] S. de Groot, P., Mazur: Non equilibrium thermodynamics, Dover, 1984, (reprint). 

[DLS01] B. Derrida, J.L. Lebowitz, E. Speer: Free Energy Functional for Nonequilibrium Systems: 
An Exactly Solvable Case, Physical Review Letters, 2001. 

[ECM93] D.J. Evans, E.G.D. Cohen, G., Morriss: Probability of second law violations in shearing 
steady flows, Physical Review Letters, 70, 2401-2404, 1993. 



2/ febbraio/2008; 17:31 



15 



[EM90] D.J. Evans and G.P. Morriss, Statistical Mechanics of Nonequilibrium Fluids, Academic 
Press, 1990, New- York. 

[EPR99] J. P. Eckmann, C.A. Pillet, L. Rey Bellet: Non-Equilibrium Statistical Mechanics of An- 
harmonic Chains Coupled to Two Heat Baths at Different Temperatures, Communications in Math- 
ematical Physics, 201, 657 - 697, 1999. 

[ES93] D.J. Evans, S. Sarman: Equivalence of thermostatted nonlinear responses, Physical Review, 
E 48, 65-70, 1993. 

[EY05] J. P. Eckmann, L.S. Young: Nonequilibrium energy profiles for a class of ID models, Com- 
munications in Mathematical Physics, ??, ????, 2005. 

[Ga96] G. Gallavotti: Chaotic hypothesis: Onsager reciprocity and fluctuation-dissipation theorem, 
Journal of Statistical Physics, 84, 899-926, 1996. 

[Ga98] G. Gallavotti, Chaotic dynamics, fluctuations, non- equilibrium ensembles, Chaos, 8, 384- 
392, 1998. 

[Ga99] G. Gallavotti: Statistical Mechanics, Springer Verlag, Berlin, 1999. 

[Ga04] G. Gallavotti, Entropy production in nonequilibrium stationary states: a point of view, 
Chaos, 14, 680-690, 2004. 

[GBG04] G. Gallavotti, F. Bonetto, G. Gentile: Aspects of the ergodic, qualitative and statistical 
theory of motion, Springer- Verlag, p. 1-434, 2004. 

[GC95] G. Gallavotti, E.G.D. Cohen: Dynamical ensembles in non- equilibrium statistical mechanics, 
Physical Review Letters, 74, 2694-2697, 1995. 

[GR97] G. Gallavotti, D. Ruelle: SRB states and nonequilibrium statistical mechanics close to 
equilibrium, Communications in Mathematical Physics, 190 (1997), 279-285. 

[GL03] S. Goldstein, J. Lebowitz, On the (Boltzmann) entropy of Nonequilibrium systems, Physica 
D, 193, 53-66, 2004. 

[Ku98] J. Kurchan, Fluctuation theorem for stochastic dynamics, Journal of Physics A, 31, 3719- 
3729, 1998. 

[Le93] J.L. Lebowitz: Boltzmann's entropy and time's arrow, Physics Today, 32-38, 1993 (Septem- 
ber). 

[LS99] J.L. Lebowitz, H. Spohn: The Gallavotti- Cohen Fluctuation Theorem for Stochastic Dy- 
namics, Journal of Statistical Physics, 95, 333-365, 1999. 

[Ma99] C. Maes: The Fluctuation Theorem as a Gibbs Property, Journal of Statistical Physics, 95, 
367-392, 1999. 



2//ei)f)raio/2008; 17:31 



16 



[Ru76] D. Ruelle: A measure associated with axiom A attractors, American Journal of Mathematics 
98, 619-654, 1976. 

[Ru96] Ruelle, D.: Positivity of entropy production in nonequilibrium statistical mechanics, Journal 
of Statistical Physics, 85, 1-25, 1996. 

[Ru97a] Ruelle, D.: Smooth dynamics and new theoretical ideas in non- equilibrium statistical me- 
chanics, Journal of Statistical Physics, 95, 393-468, 1999. 

[Ru97b] D. Ruelle: Entropy production in nonequilibrium statistical mechanics, Communications 
in Mathematical Physics, 189, 365-371, 1997. 

[Ru99a] Ruelle, D.: A remark on the equivalence of isokinetic and isoenergetic thermostats in the 
thermodynamic limit, Journal of Statistical Physics, 100, 757-763, 2000. 

[Ru99b] D. Ruelle: A remark on the equivalence of isokinetic and isoenergetic thermostats in the 
thermodynamic limit, Journal of Statistical Physics, 100, 757-763, 2000. 

[Si72] Ya.G. Sinai: Gibbs measures in ergodic theory, Russian Mathematical Surveys 166 (1972), 
21-69. 

[Si94] Ya.G. Sinai: Topics in ergodic theory, Princeton Mathematical Series, Vol. 44, Princeton 
University Press, Princeton, 1994. 



2/ febbraio/2008; 17:31 



17 



